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Abstract 



We study the critical behavior of the ferromagnetic Ising model on random trees as well 
as so-called locally tree-like random graphs. We pay special attention to trees and graphs 
Ph " with a power-law offspring or degree distribution whose tail behavior is characterized by its 

power-law exponent r > 2. We show that the critical temperature of the Ising model equals 
the inverse hyperbolic tangent of the inverse of the mean offspring or mean forward degree 
distribution. In particular, the inverse critical temperature equals zero when r S (2, 3] where 
this mean equals infinity. 

We further study the critical exponents S,f3 and 7, describing how the (root) magnetiza- 
£> I tion behaves close to criticality. We rigorously identify these critical exponents and show that 

they take the values as predicted by Dorogovstev, et al. [9] and Leone et al. [T7]. These values 
depend on the power-law exponent r, taking the mean-field values for t > 5, but different 
C*~) \ values for r 6 (3, 5). 

7— 1 [ 

1 Introduction 

In the past decades complex networks and their behavior have attracted much attention. In the 
^ | real world many of such networks can be found, for instance as social, information, technological 
ctf I and biological networks. An interesting property that many of them share is that they are scale 
free [22] • This means that their degree sequences obey a power law, i.e., the fraction of nodes 
that have k neighbors is proportional to k~ T for some r > 1. We therefore use power-law random 
graphs simple model for real-world networks. 

Not only the structure of these networks is interesting, also the behavior of processes living 
on these networks is a fascinating subject. Processes one can think of are opinion formation, 
the spread of information and the spread of viruses. An extensive overview of complex networks 
and processes on them is given by Newman in [22]- It is especially interesting if these processes 
undergo a so-called phase transition, i.e., a minor change in the circumstances suddenly results 
in completely different behavior. Examples of such phase transitions include the sudden East 
European revolution in 1989 [16] and the unusual swine flu outbreak in 2009 [5]. 
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Physicists have studied the behavior near phase transitions, the critical behavior, on complex 
networks for many different models, see [10] for an overview. Many of these results have not been 
mathematically rigorously proved. One of the few models for which rigorous results have been 
obtained is the contact process [I], where the predictions of physicists, in fact, turned out not to 
be correct. A mathematical treatment of other models is therefore necessary. 

We focus on the Ising model, a paradigm model for the study of phase transitions [23] EH G5] . 
In this model a spin value that can be either +1 or —1 is assigned to every vertex. These spins 
influence each other with ferromagnetic interactions, i.e., neighboring spins prefer to be aligned. 
The strength of these interactions depends on the temperature. The first rigorous study of the 
Ising model on a random graph was performed by De Sanctis and Guerra in [26], where the high 
and zero temperature regime of the Ising model on the Erdos-Renyi random graph were analyzed. 
Later, in [Bj, Dembo and Montanari analyzed the Ising model on locally tree-like random graphs 
with a finite- variance degree distribution for any temperature. In [Sj, we generalized these results 
to the case where the degree distribution has strongly finite mean, but possibly infinite variance, 
i.e., the degree distribution obeys a power-law with exponent r > 2. An analysis of the critical 
behavior, however, was still lacking. 

In this article, we rigorously study the critical behavior of the Ising model on power-law random 
graphs by computing various critical exponents. Predictions for the values of these exponents were 
given Dorogovtsev, et al. in [9] and independently by Leone et al. in [17] and we prove that these 
values are indeed correct. These exponents depend on the power-law exponent r. We prove that 
the critical exponents S, (3 and 7 take the classical mean-field values for r > 5, and hence also for 
the Erdos-Renyi random graph, but are different for r G (3,5). In [9] [T7] also the case r G (2,3) 
is studied for which the critical temperature is infinite. Hence, the critical behavior should be 
interpreted as the temperature going to infinity, which is a different problem from approaching a 
finite critical temperature and is therefore beyond the scope of this article. 

Our proofs always start by relating the magnetization of the Ising model on the random graph 
and various of its derivatives to the root magnetization of a rooted random tree, the so-called 
unimodular tree. After this, we identify the critical exponents related to the root magetization on 
the rooted random tree. As a result, all our results also apply to this setting, where only in the 
case of the regular tree, the mean- field critical exponents have been identified [2], and which we 
extend to general offspring distributions. 

2 Model definitions and results 
2.1 Ising model on finite graphs 

We start by defining Ising models on finite graphs. Consider a random graph sequence (G n )„>i. 
Here G n = (V n ,E n ), with vertex set V n = [n] = {1, . . . ,n} and with a random edge set E n . To 
each vertex i G [n] an Ising spin o~i = ±1 is assigned. A configuration of spins is denoted by 
a = (o"i)ig[„]. The Ising model on G n is then defined by the Boltzmann-Gibbs measure 
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Here, (3 > is the inverse temperature and B_ the vector of external magnetic fields B_ = (Sj) i6 [ n i G 
R n . For a uniform external field we write B instead of B_, i.e., Bi = B for all i G [n]. The partition 
function Z n ((3,B_) is the normalization constant in (12. ip . i.e., 

Z n {P,B)= exp /3 J] + (2.2) 

<T6{-1,+1}" [ (i,j)£E n ie[n] ) 

Note that the inverse temperature /3 does not multiply the external field. This turns out to be 
technically convenient and does not change the results, because we are only looking at systems at 
equilibrium, and hence this would just be a reparametrization. 

We let ( • ) denote the expectation with respect to the Ising measure a n , i.e., for every 
bounded function / : {—1, +l} n — > R, we write 

W)^= E /(*K(*)- (2.3) 

<T6{-1,+1}™ 



2.2 Thermodynamics 

We study the critical behavior of this Ising model by analyzing the following two thermodynamic 
quantities: 

Definition 2.1 (Thermodynamic quantities). For a graph sequence (G n ) n >i, 

(a) let M n ((3,B) = ^J2ie[ n ]( cr i)vn be the magnetization per vertex. Then, the thermodynamic 
limit of the magnetization per vertex equals 

M(f3,B) = lim M n (P,B). (2.4) 



m«) denote the susceptibility. Then, the ther- 
modynamic limit of the susceptibility equals 

X (P,B)= lim Xn(P,B). (2.5) 



The existence of the above limits for n — > oo has been proved in [HI Theorem 1.5], using the 
existence of the pressure per particle proved in [6] and [HI Theorem 1.4] and using monotonicity 
properties. We now define the critical temperature. We write /(0 + ) for lim^o /(^)- 

Definition 2.2 (Critical temperature). The critical temperature equals 

& = inf {/3 : M(/3, + ) > 0}. (2.6) 

Note that such a /3 C can only exist in the thermodynamic limit, but not for the magnetization 
of a finite graph, since always M n ((3, + ) = 0. The critical behavior can now be expressed in terms 
of the following critical exponents. We write f(x) x g{x) if the ratio f(x)/g(x) is bounded away 
from and infinity for the specified limit. 
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Definition 2.3 (Critical exponents). The critical exponents f3, 8, 7, 7' are defined by: 



M(0, 0+) x (p - C )P, for \ C - (2.7) 

M(f3 c , B) x B 1/s , for B\0; (2.8) 

X (/3,0 + )x(/3 c -/3H, for/3//3 c ; (2.9) 

X(P, + ) x (/3 - /3 C )"V, for \ (2.10) 



Remark. We emphasize that there is a difference between the symbol f3 for the inverse temperature 
and the bold symbol (3 for the critical exponent in (12.71) . Both uses for (5 are standard in the 
literature, so we decided to stick to this notation. 

Also note that these are stronger definitions than usual. E.g., normally the critical exponent (3 
is defined as that value such that 

M(/3,0 + ) = (/3-/3 c )^ +o ( 1 ), (2.11) 
where o(l) is a function tending to zero for (3 \ (3 C . 



2.3 Locally tree-like random graphs 

We study the critical behavior of the Ising model on graph sequences (G J n ) n >i that are assumed 
to be locally like a homogeneous random tree, to have a power-law degree distribution and to be 
uniformly sparse. We give the formal definitions of these assumptions below, but we first introduce 
some notation. 

Let the random variable D have distribution P = (pk)k>i, i-e., P[-D — k] — Pk, for k — 1, 2, 

We define its forward degree distribution by 

Hk E[D) ' V ; 

where we assume that K[D] < 00. Let K be a random variable with ¥[K = k] = and write 
v = K[K]. The random rooted tree T(D, K, £) is a branching process with I generations, where the 
root offspring is distributed as D and the vertices in each next generation have offsprings that are 
independent of the root offspring and are independent and identically distributed (i.i.d.) copies of 
the random variable K. We write T(K, t) when the offspring at the root has the same distribution 
as K. 

We write that an event A holds almost surely (a.s.) if W[A] = 1. The ball of radius r around 
vertex i, Bi(r), is defined as the graph induced by the vertices at graph distance at most r from 
vertex i. For two rooted trees T\ and T2, we write that T\ — T2, when there exists a bijective map 
from the vertices of 71 to those of T2 that preserves the adjacency relations. 

Definition 2.4 (Local convergence to homogeneous random trees). Let F n denote the law induced 
on the ball Bi(t) in G n centered at a uniformly chosen vertex i G [n]. We say that the graph 
sequence (G n ) n >i is locally tree-like with asymptotic degree distribution P when, for any rooted 
tree T with t generations 

lim ¥ n [Bi(t) ~ 71 = F[T(D,K,t) ~ T}. (2.13) 
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Note that this implies in particular that the degree of a uniformly chosen vertex of the graph 
has an asymptotic degree distributed as D. 

Definition 2.5 (Uniform sparsity). We say that the graph sequence (G n )n>i is uniformly sparse 
when, a.s., 

lim limsup - V Dil {Di> n = 0, (2.14) 

ie [n] 

where Di is the degree of vertex i and I.4 denotes the indicator of the event A. 

Note that uniform sparsity follows if ^ X^ieW ~~ ^\P\ a - s -> by the weak convergence of the 
degree of a uniform vertex. 

We pay special attention to cases where the degree distribution satisfies a power law, as defined 
in the following definition. For power-law degree distributions, not all moments of the degrees are 
finite, which has severe consequences for the critical behavior of the Ising model. 

Definition 2.6 (Power laws). We say that the distribution P = (pk)k>i obeys a power law with 
exponent r when there exist constants C v > c p > such that, for all k = 1,2,..., 

Cpk'^-V < J^Pe < Cpk-^-V. (2.15) 

£>k 

2.4 The random Bethe tree 

We next extend our definitions to the random tree T(D, K, oo), which is an infinite random tree. 
One has to be very careful in defining a Gibbs measure on this tree, since trees suffer from the fact 
that the boundaries of intrinsic (i.e., graph distance) balls in them have a size that is comparable 
to their volume. We can adapt the construction of the Ising model on the regular tree in [2] to 
this setting, as we now explain. For > 0, B 6 R, let fi^/ B be the Ising model on T(D, K, t) with 
+ respectively free boundary conditions. For a function / that only depends on T(-D, K, m) with 
m < t, we let 

(/)„+// = lim </> +// . (2.16) 

Below, we argue that these limits indeed exist and are equal for B > 0. This defines a unique 
infinite volume Gibbs measure /j^^ on the random Bethe tree. The quantity M(/3, B) is the 
expected root magnetization for this infinite volume Gibbs measure on the random Bethe tree. 
Our results also apply to this setting under the assumption that the degree of the root obeys a 
power law in (I2.15P or that E[K 3 ] < oo. The critical value C for the root magnetization is again 
defined by (1231) . 

2.5 Main results 

We now present our main results which describe the critical behavior of the Ising model on power- 
law random graphs and random trees with power-law offspring distribution. We first give an 
expression for the critical temperature: 
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Theorem 2.7 (Critical temperature). Assume that the random graph sequence (G n ) n >x is locally 
tree-like with asymptotic degree distribution P and is uniformly sparse. Then, a.s., the critical 
temperature (3 C of (G n ) n >i and of the random Bethe tree T(D, K, oo) equals 



(3 C = atanh(l/V) 



(2.17) 



Near the critical temperature the behavior of the Ising model can be described by critical 
exponents. The values of these critical exponents for different values of r are stated in the following 
theorem: 

Theorem 2.8 (Critical exponents). Assume that the random graph sequence (G n )n>i is locally 



tree-like with asymptotic degree distribution P that obeys K[K 



< oo or 



a power law with exponent 



r G (3, 5], and is uniformly sparse, or that the random Bethe tree obeys E[i^ 3 ] < oo or a power law 
with exponent r G (3,5]. Then, the critical exponents f3,d and 7 defined in Definition \2.3\ exist 
and satisfy 





re (3,5) 


E[K 3 ] < 00 




V(r -3) 


1/2 


S 


t - 2 


3 


1 


1 


1 


l' 


> 1 


> 1 



For the boundary case r = 5 there are logarithmic corrections for (3 = 1/2 and 8 = 3, but not for 
7 = 1 and for the lower bound 7' > 1. Indeed, (12.91) holds with 7 = 1 and the lower bound in 
( 12.1 Op holds with 7' = 1, while 



M(/3,0+) 



log 1/03 -/3 C 



1/2 



for 0\p c 



M(p c , B) x ( , B ,^ ) 1/3 f ( ,r B \. 0. 



log(l/5)J 



(2.18) 



Unfortunately, we cannot prove that the critical exponent 7' exists, see the discussion in the 
next section for more details on this issue. 



2.6 Discussion and open problems 

In this section, we discuss relations to the literature, possible extensions and open problems. 

The Ising model on random trees and random graphs. A key idea to analyze the Ising 
model on random graphs is to use the fact that expectations of local quantities coincide with the 
corresponding values for the Ising model on suitable random trees [6] . Statistical mechanics models 
on deterministic trees have been studied extensively in the literature (see for instance [21 [18] and 
its relation to "broadcasting on trees" in [T2J EI])- The analysis on random trees is more recent 
and has been triggered by the study of models on random graphs. Extensions beyond the Ising 
model, e.g., the Potts model, pose new challenges [7]. 
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Relation to the physics literature. Theorem 12 . 81 confirms the predictions in [9], [TT] . For r < 3, 
one has v = oo and hence fl c = by Theorem 12 .7] so that the critical behavior coincides with the 
infinite temperature limit. Since in this case there is no phase transition at finite temperature, 
we do not study the critical behavior here. For r = 5, in [9], also the logarithmic correction for 
j3 = 1/2 in (12.18p is computed, but not that of 8 = 3. 

The critical exponents 7' and other critical exponents. Theorem 12.81 only gives a lower 
bound on the critical exponent 7'. It is predicted that 7' = 1 for all r > 3, while there are also 
predictions for other critical exponents. For instance the critical exponent a.' for the specific heat 
in the low-temperature phase satisfies cx' = when E[fC 3 ] < 00 and a' = (r — 5)/(r — 3) in the 
power-law case with r G (3,5) (see [9_1 [TT]). We prove the lower bound 7' > 1 in Section I6T21 
below, and we also present a heuristic argument that 7' < 1 holds. The critical exponent a.' for 
the specific heat is beyond our current methods, partly since we are not able to relate the specific 
heat on a random graph to that on the random Bethe tree. 

Light tails. The case E[fC 3 ] < 00 includes all power-law degree distributions with r > 5, but 
also cases where P does not obey a power law. This means, e.g., that Theorem 12.81 also identifies 
the critical exponents for the Erdos-Renyi random graph where the degrees have an asymptotic 
Poisson distribution. 

Inclusion of slowly varying functions. In Definition 12. 6 \ we have assumed that the asymptotic 
degree distribution obeys a perfect power law. Alternatively, one could assume that ^2 i>k Pi x 
L(k)k^^ T ^ 1 ^ for some function k 1— > L(k) that is slowly varying at k = 00. For r > 5 and any 
slowly varying function, we still have E[i^ 3 ] < 00, so the results do not change and Theorem 12.81 
remains to hold. For r G (3,5], we expect slowly varying corrections to the critical behavior in 
Theorem 12 .81 For example, E[K 3 ] < 00 for r = 5 and L(k) = (log A;) -2 , so that the logarithmic 
corrections present for r = 5 disappear. 

Beyond the root magnetization for the random Bethe tree. We have identified the critical 
value and some critical exponents for the root magnetization on the random Bethe tree. The 
random Bethe tree is a so-called unimodular graph, which is a rooted graph that often arises as 
the local weak limit of a sequence of graphs (in this case, the random graphs {G n ) n >i). See [Tj [3] 
for more background on unimodular graphs and trees, in particular, 1~(D, K, 00) is the so-called 
unimodular Galton- Watson tree as proved by Lyons, Pemantle and Peres in [2D]. One would expect 
that the magnetization of the graph, which can be defined by 



where B${t) is the graph induced by vertices at graph distance at most t from the root and 
|i?^(t)| is the number of elements in it, also converges a.s. to a limit. However, we expect that 
Mt{(3, B) 7^ M(/3, B) due to the special role of the root 0, which vanishes in the above limit. Thus 
one would expect to believe that Mt{(3,B) equals the root magnetization of the tree where each 
vertex has degree distribution K + 1. Our results show that also Mt((3,B) has the same critical 
temperature and critical exponents as M((3,B). 




(2.19) 
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Relation to the Curie- Weiss model. Our results show that locally tree-like random graphs 
with finite fourth moment of the degree distribution are in the same universality class as the mean- 
field model on the complete graph, which is the Curie- Weiss model. We further believe that the 
Curie- Weiss model should enter as the limit of r — > oo for the r-regular random graph, in the sense 
that these have the same critical exponents (as we already know), as well as that all constants 
arising in asymptotics match up nicely (cf. the discussion at the end of Section l6\2l) . Further, our 
results show that for r G (3, 5], the Ising model has different critical exponents than the ones for 
the Curie- Weiss model, so these constitute a set of different universality classes. 

Organization of the article. The remainder of this article is organized as follows. We start with 
some preliminary computations in Section |3j In Section H] we prove that the critical temperature 
is as stated in Theorem 12.71 The proof that the exponents stated in Theorem 12.81 are indeed the 
correct values of (3 and S is given in Section 15.31 The value of 7 is identified in Section El where 
also the lower bound on 7' is proved and a heuristic is presented for the matching upper bound. 



3 Preliminaries 

An important role in our analysis is played by the distributional recursion 

h^^B + f^ah?), (3.1) 

i=l 

where 

= atanh(/3tanh(/i)), (3.2) 

with P = tanh(/3), and where = B, (K t )t>u are i-i-d. with distribution p and (hf )%>i are 
i.i.d. copies of independent of K t . In [8], Proposition 1.7], we have proved that this recursion 
has a unique fixed point h for all > and B > 0. Whenever we write h or hi this is a random 
variable distributed as the fixed point of (13. 1 1) . Since h is a fixed point, we can interchange 

h — B + ^2f =1 £,{hi) in expectations and we often do this. 

This fixed point h yields the random field acting on the root of the random Bethe tree 
T(D, K, 00) due to its offsprings. In particular we can use the fixed point h to give an explicit 
expression for the magnetization: 

Proposition 3.1 (Magnetization). Assume that the random graph sequence (G n )n>i is locally 
tree-like with asymptotic degree distribution P that obeys K[K] < 00 or a power law with exponent 
t G (2,3) and is uniformly sparse. Then, a.s., for all (3 > and B > 0, the thermodynamic limit 
of the magnetization per vertex exists and is given by 

D 

M(p,B)=E tanh(.B + ^£(/i i )) , (3.3) 

where 

(i) D has distribution P: 



i=l 
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(ii) (/ij)i>i are i.i.d. copies of the fixed point of the distributional recursion (13. lft ; 

(iii) D and {hi)i>\ are independent. 

The same holds on the random Bethe tree 7~(D, K, oo). 

This proposition was proved in J8j Corollary 1.6(a)] by differentiating the expression for the 
thermodynamic limit of the pressure per particle that was first obtained. Here we present a more 
intuitive proof: 

Proof of Proposition \3.1\ Let be a vertex picked uniformly at random from [n] and E n be the 
corresponding expectation. Then, 



M r 



(/3,5) = -V( ( T J )=E n [(^)]. (3.4) 



Denote by the expectations with respect to the Ising measure with +/free boundary con- 

ditions on vertices at graph distance i from <\>. Note that (cr^) £,+// ^ only depends on the spins of 
vertices in B^{i). By the GKS inequality [T5] , 

< fo) < (3.5) 

Taking the limit n — > oo, the ball B </,(£) has the same distribution as the random tree T(D, K, £), 
because of the locally tree-like nature of the graph sequence. Conditioned on the tree T, we can 
prune the tree, see J6J Lemma 4.1], to obtain that 

D 



(a^f = tanh(5 + ^£(/*f ■ (3.6) 

Similarly, 



i=l 



D 

,'(*+!) 



tanh(s + ^e(^ _1) )), (3.7) 



where ft,,- + also satisfies (I3.ip . but has initial value h'^ = oo. Since this recursion has a unique 
fixed point [HI Proposition 1.7], we prove the proposition by taking the limit £ — > oo and taking 
the expectation over the tree T(D, K, oo). □ 

To study the critical behavior we investigate the function = atanh(/3tanhx) and prove 
two important bounds that play a crucial role throughout this paper: 

Lemma 3.2 (Properties of x i— >■ £(x)). For all x, (3 > 0, 

fix ^—x 3 < £(x) < fix. (3.8) 

3(1 -/3 2 ) 



The upper bound holds with strict inequality ifx,f3 > 0. 



9 



x 2 



Proof. By Taylor's theorem, 

^)=£(o) + r(o)x + r(C)^r, (3-9) 

for some ( G (0,x). It is easily verified that £(0) = 0, 



/3(1 — tanh x 



1 - f3 2 tanh 2 x 



(3, (3.10) 



x=0 

and 



£ , (0 = 2/3(1 -/3 2 ) (tanh Q(l- tanh 2 Q < n 
(l-/3 2 tanh 2 C) 2 

thus proving the upper bound. If x, /3 > then also ( > and hence the above holds with strict 
inequality. 

For the lower bound, note that £"(0) = and 

(0 - J^-^-f^O h _ 3(1 - ^ tanh 2 C - ^ tanh* C 
(1 -/3 2 tanh 2 C) 3 V 

2/3(1 -/3 2 )(1- tanh 2 Q 2/3 

(1 -/3 2 ) 2 (1 -tanh 2 C) ~ 1-/3 2 

Thus, for some ( 6 (0,1 



> ~z 2 : = -- * n - (3.i2) 



e(x) = e(o) + £'(o)* + no)y + r (o^ > ^ - j^fr ( 3 -is) 

□ 

We next study tail probabilites of (pk)k>o- Here, for a probability distribution (qk)k>o on the 
integers, we write q> k = J2 e > k qe- 

Lemma 3.3 (Tail probabilities of (pk)k>o)- Assume that (12.151) holds for some r > 2. Then, for 
the size-biased distribution defined in li2.1S\) . there exist < c p < C p such that, for all k > 1, 



k- {T - 2) < p> k < Cpk-<r- 2 \ (3.14) 



Proof. The lower bound follows directly from the fact that p> k > (k + l)p>k+ifE[D], and (I2.15p . 
For the upper bound, we note that for any probability distribution (qk)k>o on the non-negative 
integers, we have the partial summation identity 

E »/(*) = /(0) + E 9>tW) - Sit - 1)], (3.15) 

fc>o e>i 

provided that either f(k)q>k ->0asfc->ooorfci-> f(k) is either non- decreasing or non-increasing. 
Indeed, 

k 

E = /(°) + E «*[/(*) - = /(°) + E » £[/o - - !)]■ ( 3 - 16 ) 

fc>0 fc>0 fc>0 £=1 
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Interchanging the summation order (which is allowed by our assumptions) provides the proof. 
We start by proving bounds on p> fc . We rewrite 

l>k [ J i>0 

where /(£) = {I + l)l {e > k} /E[D]. By (1315|) with = p m , for k > 1 so that /(0) = 0, 

p>, = Em - w - = nip + (3.18) 

From (P2H5D, it follows that 

"-^(* +1 '^ 2,+ E^^r-", (3.i 9 ) 

L J £>fc+l L 1 

so that there exists a constant C p such that 

P> k < C p hr<r- 2 \ (3.20) 

□ 

When computing the critical exponents for r G (3, 5], we often split the analysis into two cases: 
one where K is small and one where K is large. For this we need bounds on truncated moments 
of K which are the content of the next lemma. 

Lemma 3.4 (Truncated moments of K). Assume that (I2.15P holds for some r > 2. Then there 
exist constants C a ^ T > such that, as £ — » oo, 

r a i (C ar £ a -( T -V whena>r-2, 

E \K a l iK<£} ] < { a T 3.21 

L {Ksm - \ Cr 2rlogi w hena = r-2. V ; 

and, when a < r — 2, 

E [K a t {K>e} ] < C a , r r~ (T ~ 2) . (3.22) 
Finally, when r = 5, there exists a constant > such that, as I — >■ oo, 

E [K{K -1){K- 2)t {Km ] > c 3 , 5 log£. (3.23) 

Proof. We start by bounding the truncated moments of K. We rewrite, using (13.151) and with 
f(k) = k a l {k < e} , 

oo oo L^J 

e [K a i {K < e} ] = = £)[/(*) - f( k - l )\p>* < X> a - - l T\p&- ( 3 - 24 ) 

k=0 k=l k=l 

Using k a — (k — l) a = a A, x a ~ l dx < a/c a_1 , we arrive at 

E [K a l {K < e} ] < aC p k a - x k- {T -^ < aC p ^ k a ~( T ~ l \ (3.25) 

k=l k=l 
11 



Note that k y-> k a ^ 1 - ) is either increasing or decreasing. Hence, 

W+ 1 f £+2 

k a ~ {T - l) < I k a - {T -^dk. (3.26) 

fe=i 



whereas for a = r — 2, 



For a > t - 2, 

e+2 2 

~ a + 2-r v ; 

pl+2 

J k a - { - T - 1) &k<2\ogl. (3.28) 

Similarly, for a < t — 2, 

E [iT a ll { x>n] = m>>^ + $> a - - l)>>fc (3.29) 

oo 

< C p [£] a - (r - 2) + aC p ka '\ k + l)" (r " 2) < C a ,X _(T_2) . 
Finally, we prove (I3.23p . for which we compute with f(k) = k(k — l)(k — 2), 

oo £ oo i 

E [K(K -1)(K- 2)11 {X < £} ] = - f( k ~ !)1 E Pi = E 3 (* "!)(*- 2 ) E ^ ( 3 ' 3 °) 

fc=l i=fc fc=3 «=fc 

We bound this from below by 

E [K{K -1)(K- 2)l {K < e} ] > ^ 3(A; - - 2 )[P>fc " P>*1- ( 3 - 31 ) 

fc=0 



By Lemma [3.31 for r = 5, the contribution due to p>i is at most 

i 3/2 p>£ < C p t~^ 2 = o(l), (3.32) 
while the contribution due to p>k and using 3(k — l)(k — 2) > k 2 for every k > 4, is at least 

c P J2 k ^> c P / — = c p [log(V7+l)-log4], (3.33) 
which proves the claim by chosing the constant 03,5 correctly. □ 

4 Critical temperature 

In this section we compute the critical temperature. 



Proof of Theorem pT7[ Let /3* = atanh(l/z/). We first show that if f3 < /?*, then 



KmM(f3,B) = 0, (4.1) 

which implies that (3 C > (3*. Later, we show that if \im.B\oM(/3,B) = then /3 < /?*, implying 
that & < /?*. 
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Proof of fl c > fl*. Suppose that fl < fl* . Then, by the fact that tanhx < x and Wald's identity, 



M{fl,B) = E 



tanh \B + ^£(hi)j 



<B + E[D]E[£(h)]. (4.2) 



We use the upper bound in Lemma 13.21 to get 

E[£(/i)] = E[atanh(/3tanh/i)] < flE[h] = fl (B + vE[£(h)]) . (4.3) 

Further, note that 

E[£(h)] = E[atanh(/3tanh/i)] < fl, (4.4) 

because tanh/i < 1. Applying inequality (14. 3p £ times to (14. 2p and subsequently using inequal- 
ity (14. 4p once gives 

1 _ (RuY 

M(P, B)<B + Bf3E[D] K-^- + (3E[D] (f3uf ■ (4.5) 

1 — flu 



Hence, 



/ - 1 - (BuY 
M((3, B) < limsup B + Bf3E[D] + f3E[D](f3u) 1 

t-^OQ \ 1 — flU 

= B(l + PE[D]-X r ), (4.6) 



I- flu 

because fl < fl* = Therefore, 

to M(fi, B) < hm B (l + emjZfi) = 0- (4-7) 
This proves the lower bound on fl c . 

Proof of fl c < fl* . We adapt Lyons' proof in [18] for the critical temperature of deterministic 
trees to the random tree to show that fl c < fl*. Assume that lims\o M(fl, B) = 0. Note that 
Proposition 13 . II shows that the magnetization M(fl, B) is equal to the expectation over the random 
tree T(D, K, oo) of the root magnetization. Hence, if we denote the root of the tree T(-D, K, oo) 
by <p, then it follows from our assumption on M(fl,B) that, a.s., lims\ j o( cr </>) — 0. 

We therefore condition on the tree T = T(D, K, oo) and if we suppose that lim B \ jO ( cr 0) = 0, 
then also lim^o h((f)) = 0. Because of (13. ip . we must then have, for all v G T, 

lim lim h e ' + (v ) = 0, (4.8) 

where we can take + boundary conditions, since the recursion converges to a unique fixed point [H 
Proposition 1.7]. Now, fix < flo < fl and choose £ large enough and B small enough such that, 
for some e = e(flo, fl) > that we choose later, 

h e ' + (v) < s, (4.9) 
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for all v G T with \v\ = 1, where \v\ denotes the graph distance from <fr to v. Note that h e,+ (v) = 
oo > e for v G T with |t> | = £. 

As in [18] . we say that II is a cutset if II is a finite subset of T \ {<p} and every path from <fi 
to infinity intersects II at exactly one vertex v G II. We write v < II if every infinite path from v 
intersects II and write a < IT if a < IT and a ^ II. Furthermore, we say that w <— v if {w, v} is an 
edge in T and \w\ = \v \ + 1. Then, since h ,+ (v) = oo > e for v G T with |u| = £, there is a unique 
cutset II, such that h e > + (v) < e for all t> < II, and for all v G II there is at least one w «— v such 
that h l > + {w) > e. 

It follows from the lower bound in Lemma [3.21 that, for v < II, 

while, for v G II, 

h £ ' + (v) = B+J2 C(tanh/i(w)) > £(e). (4.11) 
If we now choose e > such that 

2 

/3fl ^r—)=$Q, (4.12) 

V 3(1 - /3 2 ) 2 ^ 

which is possible because (3 < (3, then, 

^ + (0)>^/?i^(e). (4.13) 

Since £(e) > and lims\o lim^oo h e,+ (cj)) = 0, 

inf^^=0. (4.14) 

«en 

From [IH1 Proposition 6.4] it follows that j3 < This holds for all (3 < (3, so 

p < atanh(l/i/) = (3*. (4.15) 

This proves the upper bound on (3 C , thus concluding the proof. □ 

We next show that the phase transition at this critical temperature is continuous: 
Lemma 4.1 (Continuous phase transition). It holds that 

lim R[£(h(Pc, B))} = 0, and lim E[£(h(fi, 0+))] = 0. (4.16) 

Proof. Note that lime\olE[^(/i(/3 c , B))] = c exists, because B i-> E[£(/i(/5 c , B))] is non-decreasing 
and non-negative. Assume, by contradiction, that c > 0. By the recursion in (13. ip . for B > 0, 

K 

E[£(/i(/3, £))] = E[e(s + ^(/i,(/3, 5)))] < e (B + z/E[£(/i(/3, 5))]) , (4.17) 

8=1 
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where the inequality holds because of Jensen's inequality and the concavity of h t- > C,(h). Hence, 
c=hmE{Z(h(f3c,B))] < \xmt(B + vE[£(h(p e ,B))])=Z(vc). (4.18) 

Since £ (x) < (3 c x for x > by Lemma 13.21 and using (3 C = 1/u, we obtain 

£(i/c) < (3 c uc = c, (4.19) 

leading to a contradiction. 

An adaptation of this argument shows the second statement of the lemma. Again (3 y 
E[£(/i(/3, + ))] is non-decreasing and non-negative and we assume that linigx^ E[£(/i(/3, + ))] = 
c > 0. Then, 

c = to ERMJ, 0+))] = to to E[{(B + f){(W,B))); 
<limjme(5 + ,E[eW5))])=CM, (4.20) 
leading again to a contradiction when c > 0. □ 



5 Critical exponents: Magnetization 

In this section we prove that the critical exponents related to the magnetization, i.e., (3 and S, 
take the values stated in Theorem 12.81 The analysis involves Taylor expansions performed up to 
the right order. By these Taylor expansions, higher moments of £(h) appear. Therefore, we first 
bound these higher moments of £(h) in terms of its first moment in Section I5TT1 

In Section [531 we use these bounds to give appropriate bounds on E[£(/i)] which finally allow 
us to compute the critical exponents (3 and S in Section 15.31 



5.1 Bounds on higher moments of 

Throughout Section [5] we assume that B is sufficiently close to zero and (3 C < (3 < (3 C + e for e 
sufficiently small. We write Cj, C^i > 1 for constants that only depend on /3 and moments of K, 
and satisfy 

< liminf dip) < limsupC;(/3) < oo. (5.1) 

Here Cj appears in upper bounds, while Cj appears in lower bounds. Furthermore, we write e^, i > 1 
for error functions that only depend on /3, B, E[£(/i)] and moments of K, and satisfy 

lim sup ei(/3, B) < oo and lim ej(/3 c , B) = 0. (5.2) 

Finally, we write = K[K(K — 1) • • • (K — k + 1)] for the kth factorial moment of K, so that 
V\ = v. 
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Lemma 5.1 (Bounds on second moment of £(h)). Let (3 > (3 C and B > 0. Then, 

' C 2 E[£(h)} 2 + Be 2 when E[K 2 } < oo, 



C 2 E[£{h)] 2 log (l/E[£(/i)]) + Be 2 when r = 4, 
k C 2 E[£(h)} T - 2 + Be 2 when r e (3, 4). 



(5.3) 



Proof. We first treat the case E[K 2 ] < oo. We use Lemma I3~2l and the recursion in (13. ip to obtain 



K 



E[£(/i) 2 ] < /3 2 E[/i 2 ] = /3 2 e[(£? + ^(hi) 

i=i 

= /3 2 (B 2 + 2BvE[i{h)\ + y 2 E[^{h)} 2 + uE[^(h) 2 ]) . (5.4) 
Since 1 — /3 2 z/ > 0, because j3 is sufficiently close to (3 C and /3 C = \jv < 1, the lemma holds with 



Co 



^2 



and e 2 



5/3 2 + 2/3 2 z/E[£(/i)] 



(5.5) 



1 - /3 2 is " 1 - (5 2 v 

It is not hard to see that ( 15. 1ft holds. For e 2 the first property of ( 15.2ft can also easily be seen. The 
second property in (15. 2p follows from Lemma 14.11 

If t < 4, then EffC 2 ] = oo and the above does not work. To analyze this case, we apply the 
recursion (13. ip and split the expectation over K in small and large degrees: 

K K 



EUhf 



E 



HB + J2ahi)) 1 



-{K<i} 



i=l 



E 



-{K>e} 



i=l 



We use Lemma 13.21 to bound the first term as follows: 

K K 

?2 



E 



^(B + J2^hi)) 1 



-{K<£} 



< (3 Z E 



B + J2^ h i)) 1 



1 {K<1} 



(5.6) 



(5.7) 



i=l 



8=1 



< (3 2 (B 2 + 2BuBE[i(h)} + E[K 2 l {K < £} }E[£(h)} 2 + uE[^(hf 



For r G (3,4), 



E[K 2 t {K < e} ] < C 2 , 



«4-r 



by Lemma [3.41 while for r = 4, 



E{KH {K < 1} ] <C 2A \og> 



To bound the second sum in (15. 6p . note that £(h) < (3. Hence, 



K 



eU(b + J2^)) 1 {K>1} < P 2 E[t {K>l} } < c , T p 



i=l 



(5.8) 
(5.9) 

(5.10) 



The optimal bound (up to a constant) can be achieved by choosing i such that £ 4-T E[£(/i)] 2 and 
£ 2 ~ T are of the some order of magnitude. Hence, we choose I = l/E[£(/i)]. Combining the two 
upper bounds then gives the desired result with 

1 



Co 



l-(3 2 u 



C 2 , T {3 2 + C 0iT f3< 



(5.11) 
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where we have also used that E[£(/i)] 2 < E[£(h)} 2 log(l/E[£(/i)]), and 

B/3 2 + 2(3 2 vE[£{h)] 



e-2 



We next derive upper bounds on the third moment of £,{h): 
Lemma 5.2 (Bounds on third moment of £,{h)). Let (3 > (3 C and B > 0. Then, 

' C 3 E[£(h)} 3 + Be 3 when E[K 3 } < oo, 



mh) 3 } < { 



C 3 E[^(h)} 3 log (l/E[£(/i)]) + 5e 3 when r = 5, 



k C 3 E[£(h)} T - 2 + 5e 3 when r G (3, 5). 

Proof. For E[i^ 3 ] < oo we bound, in a similar way as in Lemma [5.11 

E[£(/i) 3 ] < /3 3 ( B 3 + 3B 2 vE[£(h)] + Wv 2 E[i{h)} 2 + 3£z/E[£(/i) 5 



(5.12) 



□ 



(5.13) 



(5.14) 



v 3 E[i{h)f + 3z/ 2 E[£(/i)]E[£(/i) 2 ] + ^vE[i{hf 



Using (15. 3p . we indeed get the bound 

E\i{hf] <C 3 E[ah)? + Be 3 , 

where 

/3 3 



l-/3 3 z/ 



(v 3 + 3i/ 2 C 2 ) 



and 



e 3 



{B 2 + 35z/e 2 + 3 {Bp + z/ 2 e 2 ) E[£(/i)] + 3 (u 2 + uC 2 ) E[£(h)} 2 } . 



(5.15) 
(5.16) 

(5.17) 



1 -f3 3 v 

For r G (3,5], we use the recursion (13.11) . make the expectation over K explicit and split in 
small and large values of K to obtain 

K K 



E[^) 3 ]=E(B + ^$,) 1 



Htf<Li/*KWLI} 



i=l 



+ E 



HB + J2tthi)) 1 



l{AT>Ll/EK(h)]J} 



(5.18) 



We bound the first sum from above by 

k 



/3 3 E 



i=l 



:<Li/EK(ft)]J} 



/3 3 (5 3 + 35 2 E[X%< Ll/EK(/l)]J} ]E[£(/i)] + 35E[^(ir - l)ll {x < u/EK(/l)]J} ]E[£(/i) 2 ] 
+ 3BE[ia {i ,< Ll/EKW]J} ]E[£(/i) 2 ] + E[tf(tf - 1)(K - 2)t {K ^ l/mhm }E[^h)} 3 
+ 3E[K(K - l)^K< im mi}Mt( h )M( h )*] + E[iO {K < LVEK(/l)]J} ]E[£(/i) 3 ] Y 
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E[K 3 t {K ^ [1/mh)]n ] < C 3 , T E[£(h)] 



(5.19) 



while, for r = 5, 

E[K 3 l {K ^ [1/mm} ] < C 3i5 (1 + log (l/E[£(/i)])) . (5.20) 
Similarly, by Lemma [3.41 for r G (3,4), 

E[K 2 l {K < ll/mhm ] < C 2)T E[^h)} T -\ (5.21) 

while, for r = 4, 

E[K 2 l {K < [1/mhm ] < C 2A (1 + log (1/E[f (fc)])) . (5.22) 

For the other terms we can replace the upper bound in the sum by infinity and use the upper 
bound on E[£(/i) 2 ] of Lemma IBTTl For the second sum in (15.181) we bound £(x) < 0, so that this sum 
is bounded from above by Co ir E[£(/i)] T ~ 2 . Combining these bounds gives the desired result. □ 



5.2 Bounds on first moment of £(h) 

Proposition 5.3 (Upper bound on first moment of £{h)). Let (3 > (3 C and B > 0. Then, there 
exists a C\ > such that 



where 



For t = 5, 



E[f (/i)] < PB + j3uE[Z(h)] - dE[^h)] s , 

3 when E[K 3 } < oo, 
t — 2 when r G (3, 5]. 



mw] <pb+ puE[t(h)} - cmm 3 log (i/E[e(/i)i) . 

Proof. We first use recursion (13.11) and rewrite it as 

K 

E[£(h)]=Ek(B + ^Z(hij)] =PB + $uE[£(h)]+T 1 + T 2 



where 



and 



Ti = E 



To = E 



£ (S + KE[£(/0] - /3 (5 + KE[Z(h)]) 



K 



t (B + Y,t(hi) )-t{B + KR[t{h)]) 



i=l 



(5.23) 
(5.24) 
(5.25) 

(5.26) 

(5.27) 
(5.28) 



Here, Ti can be seen as the error of a first order Taylor series approximation of £ (B + !OL[£(/i)]) 
around 0, whereas T 2 is the error made by replacing £(/ij) by its expected value in the sum. By 
concavity of x i— > £(x), both random variables in the expectation of 7\ and T 2 are non-positive. 
In particular, T 2 < 0, which is enough for our purposes. We next bound 7\ in the cases where 
E[i^ 3 ] < oo, r G (3, 5) and r = 5 separately. 
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Bound on T\ when E[/i 3 ] < oo. To bound Ti for E[fC 3 ] < oo we use that, a.s., 

i (B + KE[£(h))) -HB + m_i{h)}) < 0, (5.29) 
which follows from Lemma 13.21 Hence, 
Ti < E 



£ (B + KE[t(h)]) - (5 (B + KE[t(h)]) ) 1 {B+Kmh)] <^ nh i } \ ■ (5.30) 
Since £"(0) = 0, it follows from Taylor's theorem that, a.s., 

£ (B + KE[£(h))) -0(B + KE[£(h)]) = (B + KE^k)]) 3 , (5.31) 

o 



for some ( G (0, B + KE[£(h)]). UB + KE[£(h)] < atanhf, then 

r (C) = _ 2 ^ 1 -^)( 1 ~ tanh2 f 1 - 3(1 - /3 2 ) tanh 2 C - /3 2 tanh 4 () < --/3(1 - /3 2 ). (5.32) 
K) (l-/3 2 tanh 2 C) 3 ^ ^" 32^ 1 ; 

Hence, 

Ti < —0(1 ~ [(B + ^E[£(/>)]) 3 l {B+mm] < atanh i } _ 

< - A^i _ /3 2 )E[^ 3 l^ EK(/l)] < atanh i_ B} ]E[e(/ i )] 3 . (5.33) 

Bound on Ti when r G (3, 5]. For r G (3, 5], we make the expectation over K explicit: 

oo 

T i = J2 p k ( B + fcE KW]) -Hb+ , (5-34) 

fc=0 

where it should be noted that all terms in this sum are negative because of Lemma 13.21 Define 
f(k) = £(B + m[£(h)]) -$(B + kE[£(h)]) and note that f(k) is non- increasing. We use (13151) 
and Lemma [3.31 to rewrite 

oo 

Tx = = f(0)+Y,[f(k)-f(k-l)] P > k < f(0)+c p J2[f(k)-f(k-l)](k+l)^- 2 \ (5.35) 

k=0 k>l k>l 

Then, use (13 . 15[) in reverse to rewrite this as 

Tx < f(0)+c p J2f(k)[k-^ 2) -(k+l)-^] < f^){l-c P Y,k^)Hr-l)c P Y,f(k){k+l)- {T - l) . 

k>0 k>l k>0 

(5.36) 

Hence, with e = /(0)(1 - c p J2 k >i k- T )/B, 

oo 

Tx < eB + (r - l)c p (E^/i)])^ 1 + lJE^/i)])-^ (5 + fcE[£(/i)]) - /? (S + kE[Z(h)])] 

b/m(h)] 

< eB+ ( T - i) Cp mahw 1 E (« E KW]) _( ^ 1) ( B + - (b + miah)))) , 

fc=a/E[£(/i)] 

(5.37) 
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where we choose a and b such that < a < b < oo. We use dominated convergence on the above 
sum. The summands are uniformly bounded, and E[£(/i)] — > for both limits of interest. Further, 

when fcE[£(/i)] = y, the summand converges pointwise to y^^ T ^ (B + y) — (3 {B + y)J . Hence, 

we can write the sum above as 

where o(l) is a function tending to zero for both limits of interest [TJJ 216 A]. The integrand is 
uniformly bounded for y G [a, b] and hence we can bound the integral from above by a (negative) 
constant —I for B sufficiently small and ft sufficiently close to (3 C . Hence, 

E[£(/i)] < PB + PvE[£(h)] - (r - l)c p /E[£(/i)p 2 . (5.39) 



Logarithmic corrections in the bound for r = 5. We complete the proof by identifying the 
logarithmic correction for r = 5. Since the random variable in the expectation in 7\ is non-positive, 
we can bound 



7\ < E 



£ (B + KR[£(h)]) -P(B + KE[£(h)]) l {K <e /mm 



(5.40) 



Taylor expansion h i-> £(/i) to third order, using that £(0) = £"(0) = 0, while the linear term 
cancels, leads to 

- r(c) 



Ti < E 



6 



(B + KV.[Z(h)\yt {K < e/mh)]} 



(5.41) 



for some ( G (0, On the event that K < e/E[£(/t)], we thus have that ( G (0,e), and 

£'"(0 > c e = infa. e (o )e ) C'"( x ) when e is sufficiently small. Thus, 



T x < |E [{B + ^E[£(/i)]) 3 l {Jf < E/BK(h)]} ] 

< |e[£(/i)] 3 E - - 2)t {K < e/mm ] . 



(5.42) 



When r = 5, by Lemma [3.41 E [K(K — 1)(K — 2)1{^<q.] > cs^log^, which completes the proof. 

□ 

Proposition 5.4 (Lower bound on first moment of £(/i)). Let (3 > (3 C and B > 0. Then, there 
exists a constant C 2 > such that 



where 



For t = 5, 



E[£(/i)] > /3fi + ^E[£(/i)] - Cl E[£(fc)]« - fid, 

3 when E[K 3 } < 00, 
r — 2 when r G (3, 5). 



E[£(fc)] >/3B + PvE[£(h)] - C 2 E[£(/i)] 3 log(l/E[£(/ i )]) - S ei . 
Proof. We again use the split in (15.261) and we bound Ti and T 2 . 



(5.43) 



(5.44) 



(5.45) 
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The lower bound on T\. For E[fT 3 ] < oo, we use the lower bound of Lemma [3.21 to get 

Tx > — ~ [(B + KElti(h)]) 3 } ■ (5.46) 
6{i p ) 

By expanding, this can be rewritten as 

Tt > ^-^E[K 3 ]E[^(h)] 3 - Be 4 . (5.47) 

For r G (3, 5], we first split T\ in a small K and a large part. For this, write 

h{k) =£(B + kE[£(h))) -P(B + kE[£(h)]) . (5.48) 

Then, 

T\ = EihiK)) = E [ti(K)l {K < emm ] + E [ti(K)l {K>efKm]y ] . (5.49) 
To bound the first term, we again use (15.461) : 

E ^(K)1 {K < £/EKW]} ] > - g E [(5 + KE[£(/i)]) 3 l {x<e/EK(/l)]} ] • (5-50) 

3(1 — p^) 

It is easy to see that the terms _B 3 E [1{a-> £ /e[£(/i)]}] and 3B 2 E[£(h)]E [Kt^K<e/E[^(h)]}\ that we get 
by expanding the above are of the form Be. To bound the other two terms, we use Lemma [3.41 to 
obtain, for e < 1, 



3BE[£(/i)] 2 E [i^ 2 l {x < £/EKW]} ] < < 



which are all of the form Be, and 



' WE[i(h)] 2 E [K 2 ] when r G (4, 5], 

3£C 2 ,4E[£(/i)] 2 log(l/E[£(/i)]) whenr = 4, (5.51) 
k 3BC 2 , T E[£(h)} T - 2 when r G (3,4), 



C 3 , 5 E[£(/i)] 3 log(l/E[£(/0]) when r = 5, 
E [ir 3 l { ^< £/EKW]} ] E[£(/i)] 3 < <( (5.52) 

C 3iT E[^(/i)] T - 2 when r G (3, 5). 

To bound T\ for large K, we observe that 

E [ti(K)l {K>e/Em]} ] > -$BE[l {K>emm ] - (3E[^h)]E[Kl {K>e/mh)]} ]. (5.53) 

Applying Lemma [3.41 now gives, for r G (3, 5] 

E [ti(K)l{K>e/E[((h)]}] > -^C 0>T E[e(/i)] r - 2 - /3C lir E[e(/i)] r - 2 

= -C 4 E[^(/i)] T - 2 - 5e 4 . (5.54) 
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The lower bound on T 2 . To bound T 2 , we split in a small and a large K contribution: 

T 2 = E[t 2 {K)t {K < £/mm ] + E[t 2 {K)t {K>£/mm ] = If + T>, (5.55) 

where 

t 2 (k) =nB + J2 t(hi) )-Z(B + m[£{h)]) . (5.56) 



i=l 



To bound T 2 , we note that 
so that 



t 2 (k) > -/3, 



(5.57) 



T> > -pE[t {K>e/mm ] > -C 5 E[e(/i)] (r - 2)A3 , (5.58) 
where we have used Lemma [3741 in the last inequality and the Markov inequality when E[i^ 3 ] < 00. 



To bound T 2 -, we start from 



Tf = E 



£"(C) 



A' 



^(h^-KE^h)]) 1 {K < 



e/K[€(h)]} 



i=l 



(5.59) 



We use that, for some C in between B + J2Y=i £0*) and 5 + #E[£(/i)], 

K 



no > 



2/3 



1-/3 2 



i=l 



(5.60) 



to obtain 



7f > 



/3 



A" 



A 



> 



(1-/32)2 

/3 



-E 



-{Jf< 6 /EK(ft)]} 



i=l i=l 

fiz/E [(£(/*) -E[£(/i)]) 2 ] + E[iO {A -< e/EKW]} ]E [(£(/>) -E[£(/i)]) 3 ] 



(1-^)2 

+ 2E[K 2 t {K < £/mh)]} ]E[^h)]E [(£(/>) - Eit(h)]) 2 } 

( BvEUK) 2 } + 2E[K 2 t {K < £m{h)]} }E[^h)]E[ah) 2 } + z/E[£(/i) ; 



- (1-/3 2 ) 2 
Using the bounds of Lemmas 15. II and 15.21 we get, 

' - (2E[K 2 }C 2 + C 3 u) E[£(/i)] 3 - Be 5 when E[K 3 } < 00, 



(5.61) 



(2E[K 2 }C 2 + C 3 u) E[£(/i)] 3 log(l/E[£(/i)]) - 5e 5 when r = 5, 



(5.62) 



k -tt4^ (C5 fT + C 3 u) E[ah)Y~ 2 - Be, 



when r G (3, 5), 



where C' 2t = E[K 2 ]C 2 for r G (4,5) and C' 2t = C 2 for r 6 (3,4]. Here, we have also used that 
(a) E[£(/i')] 3 < E[£(/i)] 3 log(l/E[£(/i)]) for r = 5; (b) E[£(/i)] 3 < E[£(/i)] T " 2 for r G (4,5]; and (c) 
E[i^ 2 l{A-< e /E[5(/i)]}]E[^(/?,)] < ev < v for r G (3,4]. Combining the bounds on Ti and T 2 gives the 
desired lower bound on E[£(/i)]. □ 
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5.3 Critical exponents (3 and S 

It remains to show that the bounds on E[£(/i)] give us the desired result: 
Theorem 5.5 (Values of f3 and 8). The critical exponent f3 equals 



1/2 whe„E[^]<oo, g3 

1 1/(t-3) when tG (3,5), v ; 



and the critical exponent 6 equals 



,3 when E[K 3 } < oo, 
° ~ i r-2 when r G (3,5). l °' D > 



For r = 5, 



/ R — R \ 1/2 / R \ 1/3 

M(/3 - 0+)X ( log(l/(/- ft! )) ) f0r ^ A ' (t^) f0rBN °' 

(5.65) 

Proof. We prove the upper and the lower bounds separately, starting with the upper bound. 

The upper bounds on the magnetization. We start by bounding the magnetization from 
above: 



M{fi,B) = E 



D 



<B + E[D]E[£(h)]. (5.66) 



tanh B + J2^( h i) 

V i=l 

We first perform the analysis for (3. Taking the limit B \ in (I5.23P in Proposition 15.31 yields 

E[f (M] < frntiho)] - C 1 E[e(/i )]' 5 , (5.67) 

where ho = h((3,0 + ). For /3 > f3 c , by definition, E[£(/iq)] > and thus we can divide through by 
E[£(/iq)] to obtain 

mho)] 6 - 1 < (5.68) 

By Taylor's theorem, 

(3u-l<u{l-f3l){l3-(5 c ). (5.69) 

Hence, 

/ - 2 \ l/(<5— 1) 

E[e(M<f^^J (/3-/3 c ) W - 1} - (5-70) 

Using that (3 = 1/(6- 1), 

M{fi, 0+) < E[J>] f ^ 1 "^ ! 09 -A)^, (5.71) 
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from which it easily follows that 



,. M(/3,0+) ^ 
hmsup-- —-p < oo. 



(5.72) 



We complete the analysis for (3 by analyzing r = 5. Since (15.231) also applies to r = 5, (15. 72 ft holds 
as well. We now improve upon this using (I5.25P in Proposition 15. 3[ which yields in a similar way 
as above that 

Since x !->■ l/log(l/a;) is increasing on (0,1) and E[£(/i )] < C(f3 — Pc) 1 ^ 2 for some C > 0, we 
immediately obtain that 



E[£(M] 2 < 



Pv-1 



< 



dlogCl/Efe^o)]) - CilogCl/tCGS-AO 1 /*])- 



Taking the limit of /3 \ /3 C as above then completes the proof. 

We continue with the analysis for 6. Setting (3 = (3 C in (I5.23P and rewriting gives 



Pc 



1/6 



mho)] < ( gr ) ^ 



(5.74) 



(5.75) 



with h c = h(P c ,B). Hence, 



M{p c , B)<B + E[D] 



Pc 



1/5 



(5.76) 



so that, using 1/6 < 1, 



hmsup — < oo. 



The analysis for 8 for r = 5 can be performed in an identical way as for 3. 



(5.77) 



The lower bounds on the magnetization. For the lower bound on the magnetization we use 
that 

d 2 

tanh x = —2 tanh x(l — tanh 2 x) > —2, 



dx 2 



so that 
Hence, 



tanh x > x — x . 



(5.78) 
(5.79) 



D 



M(P,B) > B + E[D)E[£(h)) -E LB + J^£(/j 



i=i 



> B + E[£>]E[£(/i)] - 5e 6 - E[D(D - l)]E[£(/i)] 2 - E[D]C 2 E[e(/i)] 2A(T_2) 

= 5 + (E[D] - e 7 )E[C(h)) - Be 6 , (5.80) 
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with a A b denoting the minimum of a and b, because E[£(/i)] converges to zero for both limits of 
interest. 

We again first perform the analysis for (3 and r / 5. We get from (I5.43P in Proposition 15.41 
that 

(\ 1/(5-1) / - \ f3 

where the last inequality holds because, by Taylor's theorem, 

Pv-I>v(l-$ 2 )(I3-I3 c ). (5.82) 

Hence, 

^ in( fm >E{D] (^)Y >0 

P\f>e (P ~ PcY \ Cl J 

For t = 5, we note that (15.451) as well as the fact that log 1/x < A £ x~ e for all x G (0, 1) and some 
A e > 0, yields that 

Then again using (I5.45P yields, for some constant c > 0, 

(\ 1/2 
.Jm ) (5 ' 85) 

once more since x H- l/(log(l/x)) is increasing. 

We continue with the analysis for S. Again, setting (3 = (3 C in (I5.43p . we get 

E[^ c )]> (5.86) 



from which it follows that 

B\0 fi 1 /^ 



Cl 

as required. The extension to r = 5 can be dealt with in an identical way as in (I5.84p - (l5.85p . This 
proves the theorem. □ 

6 Critical exponents: Susceptibility 

In this section, we study the susceptibility. In Section IBTTl we identify 7, in Section [6721 we prove a 
lower bound on 7' and add a heuristic why this is the correct value. 
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6.1 The critical exponent 7 

For the susceptibility in the subcritical phase, i.e., in the high-temperature region /3 < (3 C , we can 
not only identify the critical exponent 7, but we can also identify the constant: 

Theorem 6.1 (Critical exponent 7). For K[K] < 00 and (3 < (3 C , 

1 R[D]P 



In particular, 



and hence 



x(/3,0 + ) = - , , 
1 — up 

lim x (/3, + )(/3 c -/3) = ffl^ 
PSPc A ; 1-/32 

7 = 1. 



(6.1) 



(6.2) 



(6.3) 



Proof. The proof is divided into three steps. We first reduce the suceptibility on the random graph 
to the one on the random Bethe tree. Secondly, we rewrite the susceptibility on the tree using 
transfer matrix techniques. Finally, we use this rewrite (which applies to all (3 and B > 0) to 
prove that 7 = 1. 



Reduction to the random tree. Let <ft denote a vertex selected uniformly at random from [n] 
and let denote its expectation. Then we can write the susceptibility as 



X, 



1 n 



Note that 



Li=i 



9B j 



(6.4) 



(6.5) 



which is, by the GHS inequality [13] . decreasing in external fields at all other vertices k e [n]. 
Denote by (-)*' + /^ the Ising model with +/free boundary conditions, respectively, at all vertices at 
graph distance t from <p. Then, for all t > 1, 



Xn > Erf 



(6.6) 



By introducing boundary conditions, only vertices in the ball B^t) contribute to the sum. Hence, 
by taking the limit n — > 00 and using that the graph is locally tree-like, 



X >E 



j'eT t 



(6.7) 



where the expectation now is over the random tree T t ~ T(D, K, t) with root </>. 
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For an upper bound on \ n we use a trick similar to one used in the proof of Corollary 4.5]: 
Let B'j = B if j e 5*0) and B'j = B + H j (£ B t {<j>) for some if > -5. Denote by (•)# the 
associated Ising expectation. Then, because of (I6.5p . 



{(Wi) - M( a i) 









= E 










H=0. 



By the GHS inequality, (<j^)h is a concave function of H and hence, 



E, 



9 / \ 



H=0. 



< E, 



■g {(<?<]>)h=0 - (&<!>) H=-B/2) 



Using the GKS inequality this can be bounded from above by 



E, 



H=-B/2 



E rf 



(6.8) 



(6.9) 



(6.10) 



where the equality holds because the terms depend only on the system in the ball B t (4>) and hence 
not on H . By letting n — > oo, by the locally tree-likeness, this is equal to 



^E [(fo>«.+ - (^>*./)] 



(6.11) 



where the expectation and the Ising model now is over the random tree T t ~ T(D, K, t) with root 
0. From [HI Lemma 3.1] we know that this expectation can be bounded from above by M/t for 
some constant M = M((3, B) < oo. Hence, if t — > oo, 



lim E 

t— >oo 



E 

jeT t 



< 



< lim E 

i— >oo 



E 



. (6.12) 



Rewrite of the susceptibility on trees. It remains to study the susceptibility on trees. For 
this, condition on the tree T^. Then, for some vertex j at height £ < t in the tree, denote 
the vertices on the unique path from to j by = v ,Vi, . . . ,vt = j and let, for < i < £, 
S<i = (cr Vo , . . . ,a v .). We first compute the expected value of a spin a v . on this path, conditioned 
on the spin values S<i-±. Note that under this conditioning the expected spin value only depends 
on the spin value (T Vi _ x and the effective field h Vi = hv^~^ obtained by pruning the tree at vertex 
Vi, i.e., by removing all edges at vertex V{ going away from the root and replacing the external 
magnetic field at vertex V{ by h Vi which can be exactly computed using pj Lemma 4.1]. Hence, 

{a Vi \S<W /f -- 



3<7-« < _ 1 +h« j _|_ 



-h v . 



We can write the indicators 1{<t„ =±i} = |(1 ± o'vi-i), so that the above equals 



oi 1 + ^-1 



a 



Vi-l, 



-P+h Vi _|_ Qp-fhn 



-P+h v . 



+ e- 



-P+hv, 



+ e 



1 

+ 2 



-P-h v 



-/3+h Vi 



(6.13) 



(6.14) 

*P-hv 4 \ 



+ 



/3+h Vj _|_ e -/3-h v . e -/3+h Vi _|_ e 0-h V{ 
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By pairwise combining the terms over a common denominator the above equals 

I (eP+hvi — e -P-hvi^ e -l3+h Vi _|_ e P-h v ^ _ ^ e -f3+h Vi _ e /3-h Vi ^ e /3+h Vi _|_ e ~l3-h v ^ 



a 



«*-i2 ~ (e^+K + e -^-K)(e-^ +fc «i + e /3 ~K) 

By expanding all products, this equals, after cancellations, 

e 2/3 _|_ e -2/3 e 2fc„ i _|_ e -2h Vi 



e 2/3 _|_ g-2/3 _|_ e 2ft Vi _|_ g-2h Vi e 2/3 _|_ g-2/3 + e 2ft Vi + g-2/i^ 

sinh(2/3) sinh(2/Q 
cosh(2/3) + cosh(2/i„.) cosh(2/3) + cosh(2/i„.) ' 



Using this, we have that 



/, \t,+/j _ // ic \i,+//\t,+// _ / vt,+// smh(2/3) sinh(2fe„J 

^ U^P<^-i; ; - cosh(2/3 ) + cosh(2^J + cosh(2/3) + cosh(2^J ' 

(6.17) 

Applying this recursively, we get 

^ W 1W f = {cosh(2/3) + cosh(2^J 

' / sinh(2^) " sinh(2/3) \ 

^ ^cosh(2/3) + cosh(2/i„.) fc A ^ C osh(2/3) + cosh(2fc„J / ' 

Similarly, 

sinh(2/3) 



(cr vo cr Ve ) t ' +/f = ( <7 V A a vo Yl 



J cosh(2/3) + cosh(2/i t ,J 

El sinh(2/i„.) -A- sinh(2/3) 

I cosh(2/3) + cosh(2/iJ A A cosh (2/3) + cosh (2 V) 

-q sinh(2/3) 



*,+// 



cosh (2/3) + cosh(2/i„. 



+ la Y'+lf V ( Sinh(2 ^ } FT SinH2P) I ( 6 1 9 ) 

it l v cosh(2/3) + cosh(2/i w .) ^ cos h(2/3) + cosh(2^J I ' 1 ; 



Combining the above yields 
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By taking the limit t — V oo, we obtain 



X = E 

Finally, we can rewrite 

sinh(2/3) 



sinh(2/3) 



22 i 1 < CT -o> 2 ) II coshi 2 )') : cosli(2// 1 

j^Toc i — 1 



2sinh(/3) cosh(/3) 



cosh(2/3) + cosh(2/i„.) 2 cosh(/3) 2 - 1 + cosh(2/i v .) i + cosh(2^.)-i > 

2cosh(/3) 2 



so that 



(i-K,) 2 ) E^ i n( i + 



jeToo i=i 



cosh(2/i t) J — 1 
2cosh(/3) 2 



-i 



(6.21] 



(6.22) 



(6.23) 



The rewrite in (I6.23|) is valid for all (3 and B > 0, and provides the starting point for all our results 
on the susceptibility. 



Identification of the susceptibility for j3 < (3 C . We take the limit B \ 0, for /3 < (3 C , and 
apply dominated convergence. First of all, all fields hi converge to zero by the definition of (3 C , so 



we have pointwise convergence. Secondly, 1 + 



cosh(2/i I , i ) — 1 
2cosh(/3) 2 



> 1, so that the random variable in the 



expectation is bounded from above by ^,- eToo which has finite expectation as we show below. 
Thus, by dominated convergence, the above converges to 



lim v(/9, B) = E 

B\0 



E ^ 



(6.24) 



Denote by the number of vertices at distance I from the root. Then, 



E 



E 0" 



E 



(6.25) 



because Zi > 0, a.s. Note that Z^/ (¥\D]i/ x ) is a martingale, because the offspring of the root 
has expectation K[D] and all other vertices have expected offspring v. Hence, 

hm X (/3, B) = f>[Z# = 1 + jr^y- 1 ^ = l + ?E£ (6.26) 
^=o e=i 1 ^ 

This proves (16.1 1) . We continue to prove (16. 2p . which follows by using (I5.69P and (15. 82ft : 



1 -/3 2 



f 1 - £z/ 1-/3; 



-(/3 c -/3)- 1 + l. 



(6.27) 
□ 
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6.2 Partial results for the critical exponent 7' 

For the supercritical susceptibility, we prove the following lower bound on 7': 
Proposition 6.2 (Critical exponent 7'). For r G (3,5] or E[i^ 3 ] < oo ; 

l' > 1. 



(6.28) 



Proof. We start by rewriting the susceptibility in a form that is convenient in the low-temperature 
phase. 

A rewrite of the susceptibility in terms of i.i.d. random variables. For > f3 c we start 
from (16. 23 p . We further rewrite 



X 



(/?,£) = £/^E 



(1-K) 2 ) £ exp{-£log(l 



cosh (2h Vi ) — 1 
2cosh(/3) 2 



(6.29) 



Here, and in the sequel, we use the convention that empty products, arising when i = 0, equal 1, 
while empty sums equal 0. Thus, the contribution due to £ = in the above sum equals 1. We 
write f = and Vi = oq - ■ ■ G N* for % > 1, so that the <2j_ist child of Then, 



X (/3,5) = ^ £ E 



=0 a , 



;i - Ko^iWtm ex P { - £ log ( x 

8=1 



cosh(2/i„J — 1 
2cosh(/3) 2 



(6.30) 



Let K v . be the number of children of Vi, and condition on K v . = ki for every i 6 [0,1-1], where 
we abuse notation to write [0, m] = {0, . . . , m}. As a result, we obtain that 



x (f3,B) = J2^ £ £ F(v t eT oo ,K Vi = k i Vie[0,e-i]) 

(=0 ag,...,ai_i ko,...,ki_i 



(6.31) 



x E 
Note that 



(l-K^expj-^Tlog (l 

i=l 



cosh(2/^.) — 1 
2cosh( ( 5) 2 



v t e = hVie [0,£-l] 



P(i^ = kiVi e [0,£- l],v e e To.) = Pp = fe )l{ao<feo} II P ^ = 



<fc;}- 



(6.32) 



Let Tj j be the tree that descibes all descendants of the jth child of v j, with the a^th child removed, 
and Te the offspring of v#. When v# G T^, all information of the tree T^ can be encoded in 
the collection of trees (Tij)j G [ 0i K Vi -i],ie[o,e~i] and T e , together with the sequence (a«)^~Q. Denote 
T = {(Ti,j)je[o,K Vi -i],ie[o,£-i],Ti). Then, for any collection of trees t = ((t^-^o^-i^e^-i], 

P(f = t\ K v . = h Vi G [0,£ - 1],^ G T^) = P(T = fc) JJ P(T = t itj ), (6.33) 

(jj)G[0,fci-l]x[0,£-l] 
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where the law of T is that of a Galton- Watson tree with offspring distribution K. We conclude 
that 



i-i 



X (p,B) = Y,p e E E = k )i {ao < ko} nnK = k)i {a ^ 

£=0 ao,...,a e _ 1 k ,...,k e _ 1 



(6.34) 



x E 



where (/i*(/c))f =0 satisfy the recursion relations fy£ = ft^i 



^w = 5+w +1 w) + E^m) ; 



(6.35) 



and where (^j,j)ie[o/],j>i are i-i-d. copies of the random variable h((3,B). We note that the law of 
(h*(k)Y i=0 does not depend on (aj)i 6 [o^-i], so that the summation over (aj)j 6 [ ^_i] yields 



X (/3,5) = E^ E A; P( J D = A;o)n^ P ( ir = ^) 

£=0 fco,— 



(6.36) 



i=i 



x E 



For a random variable X on the non-negative integers with E[X] > 0, we let X* be the size-biased 
distribution of X given by 

k 



n** = k) = = k). 



Then 



x(ab) = ^E^ E p(^^o)Il p ( r ^>) 



£-1 



(6.37) 



(6.38) 



=0 k ,...,k i _ 1 



i=l 



x E 



t 

(1 - ((T w ) 2 )exp { - E log ( 



1 + 



cosh(2h*(k)) - 1 
2cosh(/3) 2 



Define (/i*)l=o = {h*(D*, K*, . . . , if^, iQ)) i=0 , where the random variables (£)*, K\ , . . . , i^!, 
are independent. Then we finally arrive at 



X 



w , B) = m ± 0vfE [ (1 _ ^ exp { - i> g (l + } 

i=0 L i=l v ' 



(6.39) 
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Reduction to second moments. We now proceed towards the lower bound on 7'. Note that, 
a.s., 

{a V0 )=tanh(h*J, (6.40) 

where 

D*-l D*-l 

= b+ aKj + J2 v) < B + p + E & vo- (6-4i) 

3=1 3=1 

Therefore, 

D*-l 

(a vo ) < tanh(£? + (3 + ^ ^^j))- (6-42) 
The right hand side is independent of (h*) e i=1 , so that the expectation factorizes. Further, 

D*-l 

E[tanh(5 + p + ^ £(V?))1 ^ tanh(/3) = /3 < 1, (6.43) 

3=1 

as B \ 0, /? \ /3 C . Further, we restrict the sum over all £ to £ < m, where we take m = (f3 — /3 C ) _1 . 
This leads to 



X(P, B) > (1 ~ ^ )nD] f»'E Lp { - £ log (l + 



cosh(2/i*) - 1 
2cosh(/3) 2 



(6.44) 



We condition on all coordinates of (D*, K*, . . . , K\_ x , K e ) being at most b = ((5 — (5 c ) 3 \ which 
has probability 

F{D* <b,K*<b,..., K\_ x <b,K e <b)>{l- o(l))F(K* < b) m (6.45) 

>(l-o(l))(l-C^- 3 )) m , 

which is uniformly bounded from below by a constant for the choices m = (j3 — /3 C ) _1 and b = 
W ~ /3 c )~ 1/{t ~ 3) ■ Also, we use that f3v > 1, since (5 > (3 C . This leads us to 



X (P,B)>c x J2®b 



(6.46) 



where E fe denotes the conditional expectation given that D* < b,K* < b, . . . ,K\_ X < b,K e < b. 
Using that E[e x ] > e E W, this leads us to 



m t 1- 

X (/3,B)>c x J2^p{-J2 Eb lQ g( 1 + 

£=0 i=i 



cosh(2/i*) - 1 



}■ 



(6.47) 



2cosh(/3) 2 

Define, for a > and x > 0, the function q(x) = log ^1 + a(cosh(:r) — l)j. Differentiating leads to 

a sinh(x) 



q [x) 



1 + a(cosh(x) — 1) ' 

so that q'(x) < C q x/2 for some constant C q and all x > 0. As a result, q(x) < C q x 2 /A, so that 



(6.48) 



m 1 

Xifi, B)>c x J2 exp { - C q £I 6 [(^) 2 ] }• 

£=0 " i=l 



(6.49) 
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Second moment analysis of h*. As a result, it suffices to investigate second moments of h* 
which we proceed with now. We note that 



K*-l 



K = m + i) + B+ sV',,). (6.50) 

J'=l 

Taking expectations and using that £(ti) < /3h leads to 

E b [h*\ < P% [h* +1 ] + B + E[K* -1\K*< &]E[£(/i)]. (6.51) 

Iterating this inequality until i — i and using that E b [hf\ < B + z/E[£(/t)] (since Ef,[iT] < E[i*T]) 
leads to 



£ — % — 1 

h[h*} < p t - i {B + uE[£{h)])+ ^ p s (B + E[K* -1\K* <b]E[t;(h)}) (6.52) 



< />>-' { b + .mm + B+nK"-i\K"<mm] 

1 j3 



Similarly, 

E fo [(/i*) 2 ] < [(K +l ) 2 ] + 2PE b [h* +1 ] (B + E[K* - 1 | K* < b]E[£(h)]) (6.53) 
+ B 2 + 25E[ J ftT* - 1 | K* < b]E[£(h)] + E[(K* - 1)(K* - 2) | K* < b]E[£(h)] 2 
+ E[K* -1\K*< b}E[£(h) 2 }. 

Taking the limit B \ we thus obtain 

E 6 [(/i*) 2 ] < /3 2 E b [(^ +1 ) 2 ] + 2/31, E[iT -1\K*< b}E[£(h)} (6.54) 
+ E[(iT - l)(K* -2)\K*< b]E[£(h)] 2 + E[K* - 1 \ K* < b]E[£(h) 2 ]. 

We start analysing the case where E[i^ 3 ] < oo. By Theorem 12.81 for E[K 3 ] < oo, 

E[eW]<C (/3-/3 c ) 1/2 , (6.55) 

for some constant Cq. Substituting (I6.52p . and iterating in a similar fashion as in the proof of 
(16321 . we obtain that, for E[K 3 } < oo, 

E 6 [(^) 2 ] <C(/3-/3 c ). (6.56) 
We next extend this analysis to r e (3, 5). Note that, for every a > 0, 

E[(Kr I K* <b] = n **l 1{K ~f , (6-57) 

so that, for r G (3, 5), 

E[(iT) 2 | K* <b] < ^ (6.58) 
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Further, for r G (3, 5), 

nt(h)] < C o (0 - fa) 1 " 3 -*) , (6.59) 

and thus 

E[(K*) 2 | K* < b}E[£(h)} 2 C < b 5 - T E[£(h)} 2 < C(f3 - / g c )-(5-r)/(3-r)+2/(3-r) = _ (ggQ) 

It can readily be seen that all other contributions to E b [(/i*) 2 ] are of the same or smaller order. 
For example, when E[K 2 ] < oo and using that l/(r — 3) > 1/2 for all r G (3, 5), 

E[AT* -1\K*< b]E[£(h) 2 ] < CE[£{h)] 2 = 0{f3 - (6.61) 

while, when r G (3,4), 

E[iT -1\K*< b]E[£(h) 2 } < Cb 4 - T E[£(h)] T - 2 = C(f3 - £ c )-(4-t)/(3-t)+(t-2)/(3-t) = ^ _ 

(6.62) 

We conclude that 

I 6 [(^) 2 ] <C{fi-P c ). (6-63) 

Therefore, 

XC9, 5) > c x ]Texp { - C€(/9 - A)} = 0((/3 - /3c)" 1 ), (6.64) 

as required. 

The proof for r = 5 is similar when noting that the logarithmic corrections present in E[£(/i)] 2 
and in E[(fT*) 2 | K* < b] precisely cancel. □ 

We close this section by performing a heuristic argument to determine the upper bound on 7'. 
Unfortunately, as we will discuss in more detail following the heuristics, we are currently not able 
to turn this analysis into a rigorous proof. 

The upper bound on 7': heuristics for E[K 3 ] < 00. We can bound from above 



e=o _ i=i v ' 



(6.65) 



Now, the problem is that ftv > 1 when (5 > f} c , so that we need to extract extra decay from 
the exponential term, which is technically demanding, and requires us to know various constants 
rather precisely. Let us show this heuristically. It suffices to study large values of £, since small 
values can be bounded in a simple way. 

We blindly put the expectation in the exponential, and Taylor expand to obtain that 



We compute that 



cosh(/3) 2 = — U^. (6.67) 

1-/3 2 
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Since 

K?-l 



K^m + i + E^M, (6-68) 

3=1 

we have 

E[K* - 1], 



1 — g J E[g(fe)], (6.69) 



and 



r/t+ „ n 2/3ELRT* - ll 2 + Ef(AT* - 1)(K* -2)1(1 -5) r/1Nl9 E Hf* - 1 r/1N9l , 

E \(h* 2 ps — J L \ J A — — E £(/i) 2 + ^ — -E . 6.70 

L J (1 - /3 2 )(1 -/3) 1-/3 2 

Ignoring all error terms in the proof of Lemma 15.11 shows that 

mh)2] W f^Et^^)] 2 = «[£(/i)] 2 , (6-71) 
so in total we arrive at (also using that (3 & 1/v) 

As a result, 

E[(^) 2 ] ^u 3 (l-P)/u + 3ul/^ 



El^ti)} 2 . (6.73) 



cosh (/3) 2 l _ p 

Ignoring error terms in the computation in Lemma 15.21 shows that 

E[£(/>) 3 ] « C 3 E[e(/i)] 3 , (6.74) 

where 

C 3 = ^ (iaj + 3z/ 2 C 2 ) ps -^L- (z/ 3 + 3z/ 2 C 2 ) = ^ fz/ 3 (l - $) + 3(z/ 2 /^ 

l-/3 3 z/ 1-/3 2 (l-/3 2 )(l-/3) V 

(6.75) 

since p m 1/v. Further, again ignoring error terms in ( 15.23)) and Taylor expanding to third order 
shows that 

E[£(h)} PS PvE[£(h)] - CxEf^/i)] 3 , (6.76) 

where 

Ci = -S-±±( V C 3 + 3v 2 C 2 + u 3 ), (6.77) 
and £'"(0) = —2/3(1 — /3 2 ). Substituting the definitions for C 2 and C 3 yields 

Ct = (uC 3 + 3v 2 C 2 + v 3 ) (6.78) 

(z//3 3 z/ 3 (l - /§) + 3vP\v 2 lvf + 3v\p\\ - /3 2 ) + z/ 3 (l - /3)(1 - /3 2 )) 



3(1-/3) 

g 

3(1-/3) 



(z/ 3 (l-/3) + 3z/ 2 /3 2 ). 
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Thus, we arrive at 



mh)} 2 « (6.79) 



so that substitution into (I6.73P leads to 

M « 0, - 1) 3( r; (1 - + = 3 (j y (6 .80) 

cosh(/J) 2 ^(1/3(1-^ + 3^^) 

We conclude that 

G§i/) exp { - ^fl} < (1 + (£1/ - l))e~ 3 ^- 1 ) < e" 2 ^- 1 ). (6.81) 
cosh(p)^ J v 



This suggests that 



X (P,B) < H^Ige- 2 ^- 1 ) = O(0u — 1) _1 ), (6.82) 



v 



as required. Also, using (I6.66p . this suggests that 

lim 0u - l)x(/3, 0+) = E[ J D]/(2z/), (6.83) 

where the constant is precisely half the one for the subcritical susceptibility (see (16. ip ). It can 
be seen by an explicit computation that the same factor 1/2 is also present in the same form for 
the Curie- Weiss model. Indeed for the Boltzmann-Gibbs measure with Hamiltonian H n (a) = 
— ^J2ij< 7 i< 7 j one has (3 C = 1 and a susceptibility %(/3,0 + ) = 1/(1 — f$) for (3 < (3 C , x(/3,0 + ) = 
(1 — m 2 )/ (1 — (3(1 — m 2 )) with m the non-zero solution of m = tanh(/3m) for (3 > (3 C . Expanding 
this gives m 2 = 3(/3-l)(l + o(l)) for (3 \ 1 and hence x(P, + ) = (l + o(l))/(l-/3(l-3(/3-l))) = 
(l + o(l))/(2(/3-l)). 

It is a non-trivial task to turn the heuristic of this Section into a proof because of several 
reasons: (a) We need to be able to justify the step where we put expectations in the exponential. 
While we are dealing with random variables with small means, they are not independent, so this 
is demanding; (b) We need to know the constants very precisely, as we are using the fact that a 
positive and negative term cancel in (I6.8ip . The analysis performed in the previous sections does 
not give optimal control over these constants, so this step also requires substantial work. 

The above heuristic does not apply to r G (3,5]. However, the constant in (16.80P is always 
equal to 3, irrespective of the degree distribution. This suggests that also for r £ (3, 5], we should 
have y' < 1. 
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